Integrated multiomics implicates dysregulation of ECM and cell adhesion pathways as drivers of severe COVID-associated kidney injury.

COVID-19 has been a significant public health concern for the last four years; however, little is known about the mechanisms that lead to severe COVID-associated kidney injury. In this multicenter study, we combined quantitative deep urinary proteomics and machine learning to predict severe acute outcomes in hospitalized COVID-19 patients. Using a 10-fold cross-validated random forest algorithm, we identified a set of urinary proteins that demonstrated predictive power for both discovery and validation set with 87% and 79% accuracy, respectively. These predictive urinary biomarkers were recapitulated in non-COVID acute kidney injury revealing overlapping injury mechanisms. We further combined orthogonal multiomics datasets to understand the mechanisms that drive severe COVID-associated kidney injury. Functional overlap and network analysis of urinary proteomics, plasma proteomics and urine sediment single-cell RNA sequencing showed that extracellular matrix and autophagy-associated pathways were uniquely impacted in severe COVID-19. Differentially abundant proteins associated with these pathways exhibited high expression in cells in the juxtamedullary nephron, endothelial cells, and podocytes, indicating that these kidney cell types could be potential targets. Further, single-cell transcriptomic analysis of kidney organoids infected with SARS-CoV-2 revealed dysregulation of extracellular matrix organization in multiple nephron segments, recapitulating the clinically observed fibrotic response across multiomics datasets. Ligand-receptor interaction analysis of the podocyte and tubule organoid clusters showed significant reduction and loss of interaction between integrins and basement membrane receptors in the infected kidney organoids. Collectively, these data suggest that extracellular matrix degradation and adhesion-associated mechanisms could be a main driver of COVID-associated kidney injury and severe outcomes.

Integrative informatics analysis identifies that ginsenoside Re improves renal fibrosis through regulation of autophagy.

We previously demonstrated that ginsenoside Re (G-Re) has protective effects on acute kidney injury. However, the underlying mechanism is still unclear. In this study, we conducted a meta-analysis and pathway enrichment analysis of all published transcriptome data to identify differentially expressed genes (DEGs) and pathways of G-Re treatment. We then performed in vitro studies to measure the identified autophagy and fibrosis markers in HK2 cells. In vivo studies were conducted using ureteric obstruction (UUO) and aristolochic acid nephropathy (AAN) models to evaluate the effects of G-Re on autophagy and kidney fibrosis. Our informatics analysis identified autophagy-related pathways enriched for G-Re treatment. Treatment with G-Re in HK2 cells reduced autophagy and mRNA levels of profibrosis markers with TGF-ß stimulation. In addition, induction of autophagy with PP242 neutralized the anti-fibrotic effects of G-Re. In murine models with UUO and AAN, treatment with G-Re significantly improved renal function and reduced the upregulation of autophagy and profibrotic markers. A combination of informatics analysis and biological experiments confirmed that ginsenoside Re could improve renal fibrosis and kidney function through the regulation of autophagy. These findings provide important insights into the mechanisms of G-Re's protective effects in kidney injuries.

A large-scale retrospective study enabled deep-learning based pathological assessment of frozen procurement kidney biopsies to predict graft loss and guide organ utilization.

Lesion scores on procurement donor biopsies are commonly used to guide organ utilization for deceased-donor kidneys. However, frozen sections present challenges for histological scoring, leading to inter- and intra-observer variability and inappropriate discard. Therefore, we constructed deep-learning based models to recognize kidney tissue compartments in hematoxylin & eosin-stained sections from procurement needle biopsies performed nationwide in years 2011-2020. To do this, we extracted whole-slide abnormality features from 2431 kidneys and correlated with pathologists' scores and transplant outcomes. A Kidney Donor Quality Score (KDQS) was derived and used in combination with recipient demographic and peri-transplant characteristics to predict graft loss or assist organ utilization. The performance on wedge biopsies was additionally evaluated. Our model identified 96% and 91% of normal/sclerotic glomeruli respectively; 94% of arteries/arterial intimal fibrosis; 90% of tubules. Whole-slide features of Sclerotic Glomeruli (GS)%, Arterial Intimal Fibrosis (AIF)%, and Interstitial Space Abnormality (ISA)% demonstrated strong correlations with corresponding pathologists' scores of all 2431 kidneys, but had superior associations with post-transplant estimated glomerular filtration rates in 2033 and graft loss in 1560 kidneys. The combination of KDQS and other factors predicted one- and four-year graft loss in a discovery set of 520 kidneys and a validation set of 1040 kidneys. By using the composite KDQS of 398 discarded kidneys due to "biopsy findings", we suggest that if transplanted, 110 discarded kidneys could have had similar survival to that of other transplanted kidneys. Thus, our composite KDQS and survival prediction models may facilitate risk stratification and organ utilization while potentially reducing unnecessary organ discard.

Multiscale genetic architecture of donor-recipient differences reveals intronic LIMS1 mismatches associated with kidney transplant survival.

Donor-recipient (D-R) mismatches outside of human leukocyte antigens (HLAs) contribute to kidney allograft loss, but the mechanisms remain unclear, specifically for intronic mismatches. We quantified non-HLA mismatches at variant-, gene-, and genome-wide scales from single nucleotide polymorphism (SNP) data of D-Rs from 2 well-phenotyped transplant cohorts: Genomics of Chronic Allograft Rejection (GoCAR; n = 385) and Clinical Trials in Organ Transplantation-01/17 (CTOT-01/17; n = 146). Unbiased gene-level screening in GoCAR uncovered the LIMS1 locus as the top-ranked gene where D-R mismatches associated with death-censored graft loss (DCGL). A previously unreported, intronic, LIMS1 haplotype of 30 SNPs independently associated with DCGL in both cohorts. Haplotype mismatches showed a dosage effect, and minor-allele introduction to major-allele-carrying recipients showed greater hazard of DCGL. The LIMS1 haplotype and the previously reported LIMS1 SNP rs893403 are expression quantitative trait loci (eQTL) in immune cells for GCC2 (not LIMS1), which encodes a protein involved in mannose-6-phosphase receptor (M6PR) recycling. Peripheral blood and T cell transcriptome analyses associated the GCC2 gene and LIMS1 SNPs with the TGF-ß1/SMAD pathway, suggesting a regulatory effect. In vitro GCC2 modulation impacted M6PR-dependent regulation of active TGF-ß1 and downstream signaling in T cells. Together, our data link LIMS1 locus D-R mismatches to DCGL via GCC2 eQTLs that modulate TGF-ß1-dependent effects on T cells.

HCK induces macrophage activation to promote renal inflammation and fibrosis via suppression of autophagy.

Renal inflammation and fibrosis are the common pathways leading to progressive chronic kidney disease (CKD). We previously identified hematopoietic cell kinase (HCK) as upregulated in human chronic allograft injury promoting kidney fibrosis; however, the cellular source and molecular mechanisms are unclear. Here, using immunostaining and single cell sequencing data, we show that HCK expression is highly enriched in pro-inflammatory macrophages in diseased kidneys. HCK-knockout (KO) or HCK-inhibitor decreases macrophage M1-like pro-inflammatory polarization, proliferation, and migration in RAW264.7 cells and bone marrow-derived macrophages (BMDM). We identify an interaction between HCK and ATG2A and CBL, two autophagy-related proteins, inhibiting autophagy flux in macrophages. In vivo, both global or myeloid cell specific HCK-KO attenuates renal inflammation and fibrosis with reduces macrophage numbers, pro-inflammatory polarization and migration into unilateral ureteral obstruction (UUO) kidneys and unilateral ischemia reperfusion injury (IRI) models. Finally, we developed a selective boron containing HCK inhibitor which can reduce macrophage pro-inflammatory activity, proliferation, and migration in vitro, and attenuate kidney fibrosis in the UUO mice. The current study elucidates mechanisms downstream of HCK regulating macrophage activation and polarization via autophagy in CKD and identifies that selective HCK inhibitors could be potentially developed as a new therapy for renal fibrosis.

Genetic susceptibility to diabetic kidney disease is linked to promoter variants of XOR.

The lifetime risk of kidney disease in people with diabetes is 10-30%, implicating genetic predisposition in the cause of diabetic kidney disease (DKD). Here we identify an expression quantitative trait loci (QTLs) in the cis-acting regulatory region of the xanthine dehydrogenase, or xanthine oxidoreductase (Xor), a binding site for C/EBPß, to be associated with diabetes-induced podocyte loss in DKD in male mice. We examine mouse inbred strains that are susceptible (DBA/2J) and resistant (C57BL/6J) to DKD, as well as a panel of recombinant inbred BXD mice, to map QTLs. We also uncover promoter XOR orthologue variants in humans associated with high risk of DKD. We introduced the risk variant into the 5'-regulatory region of XOR in DKD-resistant mice, which resulted in increased Xor activity associated with podocyte depletion, albuminuria, oxidative stress and damage restricted to the glomerular endothelium, which increase further with type 1 diabetes, high-fat diet and ageing. Therefore, differential regulation of Xor contributes to phenotypic consequences with diabetes and ageing.

Combination CTLA4Ig and Anti-CD40 Ligand Treatment Modifies T and B Cell Metabolic Profiles and Promotes B Cell Receptor Remodeling in a Mouse Model of Systemic Lupus Erythematosus.

Systemic lupus erythematosus is a complex autoimmune disease with significant morbidity that demands further examination of tolerance-inducing treatments. Short-term treatment of lupus-prone NZB/WF1 mice with combination CTLA4Ig and anti-CD40 ligand, but not single treatment alone, suppresses disease for >6 mo via modulation of B and T cell function while maintaining immune responses to exogenous Ags. Three months after a 2-wk course of combination costimulatory blockade, we found a modest decrease in the number of activated T and B cells in both combination and single-treatment cohorts compared with untreated controls. However, only combination treatment mice showed a 50% decrease in spare respiratory capacity of splenic B and T cells. RNA sequencing and gene set enrichment analysis of germinal center (GC) B cells confirmed a reduction in the oxidative phosphorylation signature in the combination treatment cohort. This cohort also manifested increased expression of BCR-associated signaling molecules and increased phosphorylation of PLCγ in GC B cells after stimulation with anti-IgG and anti-CD40. GC B cells from combination treatment mice also displayed a signature involving remodeling of GPI-linked surface proteins. Accordingly, we found a decrease in cell surface expression of the inhibitory molecule CD24 on class-switched memory B cells from aged NZB/W mice that corrected in the combination treatment cohort. Because both a profound decrease in BCR signaling and remodeled immune cell metabolism enhance loss of tolerance in lupus-prone mice, our findings help to explain the restoration of tolerance observed after short-term combination costimulatory blockade.

LINE-1 promotes tumorigenicity and exacerbates tumor progression via stimulating metabolism reprogramming in non-small cell lung cancer.

BACKGROUND: Long Interspersed Nuclear Element-1 (LINE-1, L1) is increasingly regarded as a genetic risk for lung cancer. Transcriptionally active LINE-1 forms a L1-gene chimeric transcript (LCTs), through somatic L1 retrotransposition (LRT) or L1 antisense promoter (L1-ASP) activation, to play an oncogenic role in cancer progression. METHODS: Here, we developed Retrotransposon-gene fusion estimation program (ReFuse), to identify and quantify LCTs in RNA sequencing data from TCGA lung cancer cohort (n = 1146) and a single cell RNA sequencing dataset then further validated those LCTs in an independent cohort (n = 134). We next examined the functional roles of a cancer specific LCT (L1-FGGY) in cell proliferation and tumor progression in LUSC cell lines and mice. RESULTS: The LCT events correspond with specific metabolic processes and mitochondrial functions and was associated with genomic instability, hypomethylation, tumor stage and tumor immune microenvironment (TIME). Functional analysis of a tumor specific and frequent LCT involving FGGY (L1-FGGY) reveal that the arachidonic acid (AA) metabolic pathway was activated by the loss of FGGY through the L1-FGGY chimeric transcript to promote tumor growth, which was effectively targeted by a combined use of an anti-HIV drug (NVR) and a metabolic inhibitor (ML355). Lastly, we identified a set of transcriptomic signatures to stratify the LUSC patients with a higher risk for poor outcomes who may benefit from treatments using NVR alone or combined with an anti-metabolism drug. CONCLUSIONS: This study is the first to characterize the role of L1 in metabolic reprogramming of lung cancer and provide rationale for L1-specifc prognosis and potential for a therapeutic strategy for treating lung cancer. TRIAL REGISTRATION: Study on the mechanisms of the mobile element L1-FGGY promoting the proliferation, invasion and immune escape of lung squamous cell carcinoma through the 12-LOX/Wnt pathway, Ek2020111. Registered 27 March 2020 - Retrospectively registered.

Linking erythropoietin to Treg-dependent allograft survival through myeloid cells.

Erythropoietin (EPO) has multiple nonerythropoietic functions, including immune modulation, but EPO's effects in transplantation remain incompletely understood. We tested the mechanisms linking EPO administration to prolongation of murine heterotopic heart transplantation using WT and conditional EPO receptor-knockout (EPOR-knockout) mice as recipients. In WT controls, peritransplant administration of EPO synergized with CTLA4-Ig to prolong allograft survival (P < 0.001), reduce frequencies of donor-reactive effector CD8+ T cells in the spleen (P < 0.001) and in the graft (P < 0.05), and increase frequencies and total numbers of donor-reactive Tregs (P < 0.01 for each) versus CTLA4-Ig alone. Studies performed in conditional EPOR-knockout recipients showed that each of these differences required EPOR expression in myeloid cells but not in T cells. Analysis of mRNA isolated from spleen monocytes showed that EPO/EPOR ligation upregulated macrophage-expressed, antiinflammatory, regulatory, and pro-efferocytosis genes and downregulated selected proinflammatory genes. Taken together, the data support the conclusion that EPO promotes Treg-dependent murine cardiac allograft survival by crucially altering the phenotype and function of macrophages. Coupled with our previous documentation that EPO promotes Treg expansion in humans, the data support the need for testing the addition of EPO to costimulatory blockade-containing immunosuppression regimens in an effort to prolong human transplant survival.

Deep learning identified pathological abnormalities predictive of graft loss in kidney transplant biopsies.

Interstitial fibrosis, tubular atrophy, and inflammation are major contributors to kidney allograft failure. Here we sought an objective, quantitative pathological assessment of these lesions to improve predictive utility and constructed a deep-learning-based pipeline recognizing normal vs. abnormal kidney tissue compartments and mononuclear leukocyte infiltrates. Periodic acid- Schiff stained slides of transplant biopsies (60 training and 33 testing) were used to quantify pathological lesions specific for interstitium, tubules and mononuclear leukocyte infiltration. The pipeline was applied to the whole slide images from 789 transplant biopsies (478 baseline [pre-implantation] and 311 post-transplant 12-month protocol biopsies) in two independent cohorts (GoCAR: 404 patients, AUSCAD: 212 patients) of transplant recipients to correlate composite lesion features with graft loss. Our model accurately recognized kidney tissue compartments and mononuclear leukocytes. The digital features significantly correlated with revised Banff 2007 scores but were more sensitive to subtle pathological changes below the thresholds in the Banff scores. The Interstitial and Tubular Abnormality Score (ITAS) in baseline samples was highly predictive of one-year graft loss, while a Composite Damage Score in 12-month post-transplant protocol biopsies predicted later graft loss. ITASs and Composite Damage Scores outperformed Banff scores or clinical predictors with superior graft loss prediction accuracy. High/intermediate risk groups stratified by ITASs or Composite Damage Scores also demonstrated significantly higher incidence of estimated glomerular filtration rate decline and subsequent graft damage. Thus, our deep-learning approach accurately detected and quantified pathological lesions from baseline or post-transplant biopsies and demonstrated superior ability for prediction of post-transplant graft loss with potential application as a prevention, risk stratification or monitoring tool.

Recipient APOL1 risk alleles associate with death-censored renal allograft survival and rejection episodes.

Apolipoprotein L1 (APOL1) risk alleles in donor kidneys associate with graft loss, but whether recipient risk allele expression affects transplant outcomes is unclear. To test whether recipient APOL1 risk alleles independently correlate with transplant outcomes, we analyzed genome-wide SNP genotyping data on donors and recipients from 2 kidney transplant cohorts: Genomics of Chronic Allograft Rejection (GOCAR) and Clinical Trials in Organ Transplantation 01/17 (CTOT-01/17). We estimated genetic ancestry (quantified as the proportion of African ancestry, or pAFR) by ADMIXTURE and correlated APOL1 genotypes and pAFR with outcomes. In the GOCAR discovery set, we noted that the number of recipient APOL1 G1/G2 alleles (R-nAPOL1) associated with an increased risk of death-censored allograft loss (DCAL), independent of ancestry (HR = 2.14; P = 0.006), as well as within the subgroup of African American and Hispanic (AA/H) recipients (HR = 2.36; P = 0.003). R-nAPOL1 also associated with an increased risk of any T cell-mediated rejection (TCMR) event. These associations were validated in CTOT-01/17. Ex vivo studies of PMBCs revealed, unexpectedly, high expression levels of APOL1 in activated CD4+/CD8+ T cells and NK cells. We detected enriched immune response gene pathways in risk allele carriers compared with noncarriers on the kidney transplant waitlist and among healthy controls. Our findings demonstrate an immunomodulatory role for recipient APOL1 risk alleles associated with TCMR and DCAL. We believe this finding has broader implications for immune-mediated injury to native kidneys.

Reversible dysregulation of renal circadian rhythm in lupus nephritis.

BACKGROUND: We have found disruption of expression of major transcriptional regulators of circadian rhythm in the kidneys of several mouse models of lupus nephritis. Here we define the consequence of this disturbance with respect to circadian gene expression and renal homeostatic function in a mouse model of lupus nephritis. METHODS: Molecular profiling of kidneys from 47 young and 41 nephritic female NZB/W F1 mice was performed at 4 hourly intervals over a 24 h period. Disruption of major circadian transcriptional regulators was confirmed by qPCR. Molecular data was normalized and analyzed for rhythmicity using RAIN analysis. Serum aldosterone and glucose and urine sodium and potassium were measured at 4 hourly intervals in pre-nephritic and nephritic mice and blood pressure was measured every 4 h. Analyses were repeated after induction of complete remission of nephritis using combination cyclophosphamide and costimulatory blockade. RESULTS: We show a profound alteration of renal circadian rhythms in mice with lupus nephritis affecting multiple renal pathways. Using Cosinor analysis we identified consequent alterations of renal homeostasis and metabolism as well as blood pressure dipper status. This circadian dysregulation was partially reversed by remission induction therapy. CONCLUSIONS: Our studies indicate the role of inflammation in causing the circadian disruption and suggest that screening for loss of normal blood pressure dipping should be incorporated into LN management. The data also suggest a potential role for circadian agonists in the treatment of lupus nephritis.

Global transcriptomic changes in glomerular endothelial cells in mice with podocyte depletion and glomerulosclerosis.

Podocytes are a key component of the glomerular filtration barrier, and its dysfunction and eventual loss drive glomerular disease progression. Recent research has demonstrated the importance of podocyte cross-talk with other glomerular cells, such as glomerular endothelial cells (GECs), in both glomerular homeostasis and in disease settings. However, how GECs are affected globally by podocyte injury and loss in disease settings remains unclear. Therefore, to characterize the molecular changes occurring in GECs in response to the podocyte loss, we performed the transcriptomic profiling of isolated GECs after diphtheria toxin (DT)-mediated podocyte depletion in transgenic mice with podocyte-specific human DT receptor and endothelial-specific enhanced yellow fluorescent protein (EYFP) expression. DT administration led to nearly 40% of podocyte loss with the development of glomerulosclerosis. Differential gene expression analysis of isolated GECs in the diseased mice showed significant changes in pathways related to cell adhesion and actin cytoskeleton, proliferation, and angiogenesis, as well as apoptosis and cell death. However, quantification of EYFP + GECs indicated that there was a reduction in GECs in the diseased mice, suggesting that despite the ongoing proliferation, the concomitant injury and the activation of cell death program results in their overall net loss. The upstream regulator analysis strongly indicated the involvement of p53, TGF-ß1, and TNF-α as key mediators of the molecular changes occurring in GECs in the diseased mice. Our findings demonstrate significant molecular changes in GECs as a secondary consequence of podocyte loss and provide a valuable resource for further in-depth analysis of potential glomerular cross-talk mediators.

Identification of New Rare Variants Associated With Familial Autoimmune Thyroid Diseases by Deep Sequencing of Linked Loci.

CONTEXT: Genetic risk factors play a major role in the pathoetiology of autoimmune thyroid diseases (AITD). So far, only common risk variants have been identified in AITD susceptibility genes. Recently, rare genetic variants have emerged as important contributors to complex diseases, and we hypothesized that rare variants play a key role in the genetic susceptibility to AITD. OBJECTIVE: We aimed to identify new rare variants that are associated with familial AITD. METHODS: We performed deep sequencing of 3 previously mapped AITD-linked loci (10q, 12q, and 14q) in a dataset of 34 families in which AITD clustered (familial AITD). RESULTS: We identified 13 rare variants, located in the inositol polyphosphate multikinase (IPMK) gene, that were associated with AITD (ie, both Graves' disease [GD] and Hashimoto's thyroiditis [HT]); 2 rare variants, within the dihydrolipoamide S-succinyltransferase (DLST) and zinc-finger FYVE domain-containing protein (ZFYVE1) genes, that were associated with GD only; and 3 rare variants, within the phosphoglycerate mutase 1 pseudogene 5 (PGAM1P5), LOC105369879, and methionine aminopeptidase 2 (METAP2) genes, that were associated with HT only. CONCLUSION: Our study demonstrates that, in addition to common variants, rare variants also contribute to the genetic susceptibility to AITD. We identified new rare variants in 6 AITD susceptibility genes that predispose to familial AITD. Of these, 3 genes, IPMK, ZFYVE1, and METAP2, are mechanistically involved in immune pathways and have been previously shown to be associated with autoimmunity. These genes predispose to thyroid autoimmunity and may serve as potential therapeutic targets in the future.

DACH1 protects podocytes from experimental diabetic injury and modulates PTIP-H3K4Me3 activity.

Dachshund homolog 1 (DACH1), a key cell-fate determinant, regulates transcription by DNA sequence-specific binding. We identified diminished Dach1 expression in a large-scale screen for mutations that convert injury-resistant podocytes into injury-susceptible podocytes. In diabetic kidney disease (DKD) patients, podocyte DACH1 expression levels are diminished, a condition that strongly correlates with poor clinical outcomes. Global Dach1 KO mice manifest renal hypoplasia and die perinatally. Podocyte-specific Dach1 KO mice, however, maintain normal glomerular architecture at baseline, but rapidly exhibit podocyte injury after diabetes onset. Furthermore, podocyte-specific augmentation of DACH1 expression in mice protects from DKD. Combined RNA sequencing and in silico promoter analysis reveal conversely overlapping glomerular transcriptomic signatures between podocyte-specific Dach1 and Pax transactivation-domain interacting protein (Ptip) KO mice, with upregulated genes possessing higher-than-expected numbers of promoter Dach1-binding sites. PTIP, an essential component of the activating histone H3 lysine 4 trimethylation (H3K4Me3) complex, interacts with DACH1 and is recruited by DACH1 to its promoter-binding sites. DACH1-PTIP recruitment represses transcription and reduces promoter H3K4Me3 levels. DACH1 knockdown in podocytes combined with hyperglycemia triggers target gene upregulation and increases promoter H3K4Me3. These findings reveal that in DKD, diminished DACH1 expression enhances podocyte injury vulnerability via epigenetic derepression of its target genes.

Molecular Analysis of the Kidney From a Patient With COVID-19-Associated Collapsing Glomerulopathy.

Recent case reports suggest that coronavirus disease 2019 (COVID-19) is associated with collapsing glomerulopathy in African Americans with apolipoprotein L1 gene (APOL1) risk alleles; however, it is unclear whether disease pathogenesis is similar to HIV-associated nephropathy. RNA sequencing analysis of a kidney biopsy specimen from a patient with COVID-19-associated collapsing glomerulopathy and APOL1 risk alleles (G1/G1) revealed similar levels of APOL1 and angiotensin-converting enzyme 2 (ACE2) messenger RNA transcripts as compared with 12 control kidney samples downloaded from the GTEx (Genotype-Tissue Expression) Portal. Whole-genome sequencing of the COVID-19-associated collapsing glomerulopathy kidney sample identified 4 indel gene variants, 3 of which are of unknown significance with respect to chronic kidney disease and/or focal segmental glomerulosclerosis. Molecular profiling of the kidney demonstrated activation of COVID-19-associated cell injury pathways such as inflammation and coagulation. Evidence for direct severe acute respiratory syndrome coronavirus 2 infection of kidney cells was lacking, which is consistent with the findings of several recent studies. Interestingly, immunostaining of kidney biopsy sections revealed increased expression of phospho-STAT3 (signal transducer and activator of transcription 3) in both COVID-19-associated collapsing glomerulopathy and HIV-associated nephropathy as compared with control kidney tissue. Importantly, interleukin 6-induced activation of STAT3 may be a targetable mechanism driving COVID-19-associated acute kidney injury.

Key driver genes as potential therapeutic targets in renal allograft rejection.

Acute rejection (AR) in renal transplantation is an established risk factor for reduced allograft survival. Molecules with regulatory control among immune pathways of AR that are inadequately suppressed, despite standard-of-care immunosuppression, could serve as important targets for therapeutic manipulation to prevent rejection. Here, an integrative, network-based computational strategy incorporating gene expression and genotype data of human renal allograft biopsy tissue was applied, to identify the master regulators - the key driver genes (KDGs) - within dysregulated AR pathways. A 982-meta-gene signature with differential expression in AR versus non-AR was identified from a meta-analysis of microarray data from 735 human kidney allograft biopsy samples across 7 data sets. Fourteen KDGs were derived from this signature. Interrogation of 2 publicly available databases identified compounds with predicted efficacy against individual KDGs or a key driver-based gene set, respectively, which could be repurposed for AR prevention. Minocycline, a tetracycline antibiotic, was chosen for experimental validation in a murine cardiac allograft model of AR. Minocycline attenuated the inflammatory profile of AR compared with controls and when coadministered with immunosuppression prolonged graft survival. This study demonstrates that a network-based strategy, using expression and genotype data to predict KDGs, assists target prioritization for therapeutics in renal allograft rejection.

Genome-wide non-HLA donor-recipient genetic differences influence renal allograft survival via early allograft fibrosis.

Donor-recipient (D-R) differences at human leukocyte antigen (HLA) loci are currently incorporated into organ sharing, allocation and immunosuppression decisions. However, while acute rejection episodes have substantially diminished, progressive histologic damage occurs in allografts and improved long-term survival remains an unrealized goal among kidney recipients. Here we tested the hypothesis that non-HLA dependent, genome-wide D-R genetic differences could contribute to unchecked alloimmunity with histologic and functional consequences, culminating in long-term allograft failure. Genome-wide single nucleotide polymorphism (SNP) array data, excluding the HLA region, was utilized from 385 transplants to study the role of D-R differences upon serial histology and allograft survival. ADMIXTURE analysis was performed to quantitatively estimate ancestry in each D-R pair and PLINK was used to estimate the proportion of genome-shared identity-by-descent (pIBD) between D-R pairs. Subsequently, quantitative measures of recipient ancestry based on non-HLA SNPs was associated with death-censored allograft survival in adjusted Cox models. In D-R pairs of similar ancestry, pIBD was significantly associated with allograft survival independent of HLA mismatches in 224 transplants. Surprisingly, pIBD and recipient ancestry were not associated with clinical or subclinical rejection at any time post-transplant. Significantly, in multivariable analysis, pIBD inversely correlated with vascular intimal fibrosis in 160 biopsies obtained less than one year which in turn was significantly associated with allograft survival. Thus, our novel data show that non-HLA D-R differences associate with early vascular intimal fibrosis and allograft survival.

Distinct peripheral blood molecular signature emerges with successful tacrolimus withdrawal in kidney transplant recipients.

Tacrolimus (Tac) is an effective anti-rejection agent in kidney transplantation, but its off-target effects make withdrawal desirable. Although studies indicate that Tac can be safely withdrawn in a subset of kidney transplant recipients, immune mechanisms that underlie successful vs unsuccessful Tac removal are unknown. We performed microarray analyses of peripheral blood mononuclear cells (PBMC) RNA from subjects enrolled in the Clinical Trials in Organ Transplantation-09 study in which we randomized stable kidney transplant recipients to Tac withdrawal or maintenance of standard immunosuppression beginning 6 months after transplant. Eight of 14 subjects attempted but failed withdrawal, while six developed stable graft function for ≥2 years on mycophenolate mofetil plus prednisone. Whereas failed withdrawal upregulated immune activation genes, successful Tac withdrawal was associated with a downregulatory and proapoptotic gene program enriched within T cells. Functional analyses suggested stronger donor-reactive immunity in subjects who failed withdrawal without evidence of regulatory T cell dysfunction. Together, our data from a small, but unique, patient cohort support the conclusion that successful Tac withdrawal is not simply due to absence of donor-reactive immunity but rather is associated with an active immunological process.

Publisher Correction: The HIV protease inhibitor darunavir prevents kidney injury via HIV-independent mechanisms.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.